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. We analyze the structure of heavy multitrace BPS operators in M = 1 supercon- 

' formal quiver gauge theories that arise on the worldvolume of D3-branes on an affine 

toric cone. We exhibit a geometric procedure for counting heavy mesonic operators 
^ . with given ^7(1) charges. We show that for any fixed linear combination of the ^7(1) 

^ I charges, the entropy is maximized when the charges are in certain ratios. This selects 

preferred directions in the charge space that can be determined with the help of a piece 
of string. We show that almost all heavy mesonic operators of fixed U (1) charges share 
a universal structure. This universality reflects the properties of the dual extremal 
black holes whose microstates they create. We also interpret our results in terms of 
typical configurations of dual giant gravitons in AdS space. 
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1 Introduction 



Superconformal field theories (SCFTs) living on the world- volume of D3-branes transverse to 
the conical singularity of an afiine cone over a Sasaki- Einstein manifold X are holographically 
dual to gravity on AdSs x X. In this paper we explore the structure of heavy BPS states 
in these theories, describing their universal structural properties, and relating them to dual 
extremal black holes in AdSs x X. 

In [1], generating functions for enumerating mesonic BPS operators in these SCFTs were 
derived. This counting was extended to include baryonic operators in [2[|3], but in the present 
paper we consider only zero baryon charge. The authors of [1] derived the entropy associated 
to mesonic BPS states carrying a fixed linear combination of the three U{1) charges. The 
most common and physically relevant example is states of constant R-charge R. In this 
article, we focus on toric cones and introduce a geometric construction that allows us to 
refine this counting and compute the entropy associated to any triple of charges rather than 
a fixed linear combination. By extension we are able to single out the particular triple that 
maximizes the entropy subject to any linear constraint on the charges. In particular this 
triple can be computed by finding the center of gravity of the pyramid cut of the dual toric 
cone by the constraint plane. Amusingly, this amounts to suspending the pyramid from a 
piece of string, providing a novel use for a different kind of string theory. 

We will show that most heavy mesonic operators that are BPS have a universal structure. 
For a given set of U (1) charges, we find the mean distribution of trace factors in heavy multi- 
trace operators. We then quantify the scale at which this distribution may be thought of as 
defining a typical structure for heavy operators, and spell out those features that are shared 
by almost all states. The heavy operators have gravity duals that can be interpreted in terms 
of giant gravitons [2] or in terms of dual giant gravitons [4] . This allows us to reinterpret our 
results as defining typical configurations of giant gravitons and dual giants. We also analyze 
the existence of horizons in the dual supergravity solutions. 

The organization of the paper is as follows. In Section [2] we briefly review the connection 
between toric quiver theories and toric varieties. In Section [3] we present our geometric 
construction and compute the maximal entropy vector for any fixed linear combination of 
the U{1) charges. We also show that the entropies found are not large enough to produce 
macroscopic horizons in the dual gravity description. 

In Section |U we treat the system statistically, and find the mean distribution of traces 
for the multitrace operators carrying any specified triplet of U{1) charges. We also give 
the interpretation in terms of dual giant gravitons in AdSs. Strictly speaking, since giant 
gravitons and dual giant gravitons are compact objects, their wavefunctions can spread over 
their moduli spaces. As we are analyzing these systems in terms of their typical classical 
configuration, it is important to evaluate how well the wavefunction localizes at a given 
point in the configuration space. For the special case of 1/8 BPS states in A/" = 4 Yang-Mills 
theory, this analysis is carried out in Appendix [Bl In Section [5] we analyze the dual gravity 
description, and finally close in Section [6] with a discussion. In Appendix O we point out 
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some intriguing relations to topological strings and amoebae. 

2 Review 

The AdS/CFT correspondence has been extended to a set of type IIB string backgrounds 
of the form AdSs x X, where X is a five-dimensional Sasaki-Einstein manifold [5]. The 
superconformal field theories dual to string theory on these backgrounds arise from a stack 
of D3-branes located at the conical singularity associated to the metric cone of X, and are 
quiver gauge theories. The material of the present paper applies to the most extensively 
studied class of the correspondence - when the cone over X is so-called toric. 

Recently, there has been a realization that toric geometry, quiver theories, and dimer 
models give rise to an intimate web of relations. In this section we briefly review these 
connections; for a more comprehensive review the reader is referred to [6]. With these 
prerequisites, we devise a geometric construction in the dual toric cone to analyze BPS 
states carrying specified U{1) charges. A well-known example of a toric Calabi-Yau cone 
is the complex cone over a complex surface known as the first del Pezzo surface [7j. 
Geometrically, this is the projective plane CP^ blown up at a point by a sphere. We shall 
refer to this Calabi-Yau cone as dPi and shall use this variety to illustrate this section and 
the ensuing discussions. 

2.1 Toric quiver theories and planar quivers 

A toric quiver theory consists of a set of gauge groups, U{Ni) or SU{Ni); for simplicity we 
take the groups to have equal rank: Ni = N for all i. Further, the theory contains a set of 
chiral superfields Xj, each transforming in the bifundamental representation {N,N) under 
a pair of gauge groups and trivially under the others. The gauge groups and the matter 
can be presented in a directed graph, with gauge groups as nodes and matter multiplets as 
directed arrows between the two nodes under which the multiplet transforms non-trivially. 
The dynamics of the theory is encoded in a superpotential W. In toric quiver theories the 
form of W is highly restricted: it consists of a sum of monomials, with each field appearing 
linearly, and in exactly two terms with opposite signs. In Figure [1] we portray the quiver 
diagram and the superpotential of the quiver gauge theory corresponding to dPi. 

The toric quiver data and the superpotential can be incorporated into a single graph 
called the planar quiver. This is done by embedding the quiver into a torus, and "opening 
up" the arrows between the gauge groups in such a way that each arrow in the planar quiver 
corresponds to exactly one field. This can always be done so as to ensure that each term in 
the superpotential corresponds to a face in the planar quiver, i.e., for any face the arrows 
(fields) surrounding it make up a term in the superpotential. Each face is circled by the fields 
either in a clockwise or a counter-clockwise direction; the first corresponds to a negative and 
the latter to a positive superpotential term. See Figure [1] for the planar quiver of dPi. 
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U(N) 
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U(N) 



= Tr X13XI4XI1 - Tr X.^Xi^Xi, - 
+Tr X12X23XI4XI1 — Tr X12XI3XI4XI] 
+Tr X23X34X42 — Tr XI3XI4X42 



Figure 1: Left: The quiver diagram and the superpotential corresponding to dPi. Right: 
The planar quiver 

2.1.1 Mesonic operators in quiver theories 

We are interested in multi-trace mesonic gauge invariant operators, which are of the form 



They form closed loops in the quiver and the planar quiver, with definite winding numbers 
(p, q) on the torus. The minimal mesonic operators are loops that pass through any node 
of the quiver at most once, i.e., operators that cannot be split into smaller components 
by adding traces. These generate the chiral ring of mesonic operators. Thus, any mesonic 
operator can be written as a product of these minimal loops with suitably placed traces. 

The number of minimal operators can be computed from the quiver, but not all of 
these operators are independent: we must impose the vanishing of the F-terms coming from 
the superpotential W . Therefore, the minimal operators are split into equivalence classes 
under the F-flatness conditions dW/dXi = 0. These equivalence classes are in one-to-one 
correspondence with the winding numbers on the torus, i.e., two operators are equivalent if 
and only if they share the same winding numbers {p, q) . 

Applying this to our example, from the dPi quiver of Figure [1] one can compute that 
there are 24 minimal loops, and in Table [1] we split these loops into 9 equivalence classes. 
Verifying that these are the equivalence classes under the F-term constraints is left as an 
exercise to the reader. 

2.2 Brane tilings and dimer models 

We can equivalently consider the graph dual to the planar quiver: this is obtained by re- 
placing each vertex with a face, and each arrow with a perpendicular line. This dual graph 
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# 


Minimal loops 


Winding (p, q) 


1 


V V2 vl y2 


(-25 1) 


3 


y \^l vl y2 y y2 yl yl y y2 y2 y2 
-^12-^23-^34-^415 -^12-^23-^34-^415 -^12-^23-^34-^41 


(-I5I) 


3 


y y2 y3 y2 y2 yl y y yl y2 
-^12-^23-^34-^415 -^23-^34-^42, -^13-^34-^41 


(-I5O) 


3 


V vl vl vl V vl V2 V2 V V2 V2 Vl 
-^12-^23-^34-^415 -^12-^23-^34-^415 -^12-^23^34-^41 


(0,1) 


6 


V vl V3 v2 V V2 v3 vl V2 v2 V 
A12A23A34A41, A12A23A34A41, A23A34A42 

vl vl V V vl vl V V2 V2 
-^23-^34-^425 Ai3v\34v\4]^, Jii^J^^^A.^^ 


(0,0) 


2 


V V3 V2 V3 V V2 
Al3A34A4j^, A34A42A23 


(05-1) 


1 


-^12-^23-^34-^41 


(151) 


3 


V vl V3 vl vl V2 V V V2 vl 
Al2A23A34A4;^, A23A34A42, A.i3A.^4A.4^ 


(150) 


2 


V V3 vl vl V3 V 
Al3A34A4]^, A23A34A42 


(15-1) 



Table 1: Equivalence classes of minimal loops and their windings for the dP\ quiver. 

is called a brane tiling. It has a natural bipartite structure, otherwise known as a dimer 
model [6]. If we color the nodes that correspond to negative superpotential terms in black, 
and those corresponding to positive ones in white, we see that each vertex is only connected 
to vertices of opposite color. See the brane tiling of dPi as an explicit example in Figure [21 
note that we have marked the two nontrivial cycles on the torus and called them z and w. 

The advantage of working with the brane tiling is that it allows us to find the toric 
Calabi-Yau manifold on the gravity side that is dual to the quiver gauge theory. This is 
accomplished using the Kasteleyn matrix K, which is formed from the brane tiling as follows. 

First, we mark some edges in the tiling with minus-signs. This is done in such a way 
that for every face with (0 mod 4) edges surrounding it, an odd number of those faces have 
minus-signs; and for every face with (2 mod 4) surrounding edges, an even number have 
minus signs. This can always be done for any brane tiling. Then define the matrix elements 
Kij, where the index i runs over the black nodes of the tiling, and j runs over the white 
nodes. Kij is a sum over all the edges connecting the two nodes, with the following weights: 
-1, if the edge was marked with a minus sign above; z, if the edge crosses the cycle z in 
the tiling in the positive direction; z~^ if the edge crosses z in the negative direction, and 
likewise for w and w~^. The positive direction when crossing the cycles z and w can be 
chosen arbitrarily, as long as it is consistently followed for all elements. 

This is best illustrated by writing down the Kasteleyn matrix corresponding to the dPi 
brane tiling of Figure [2l In this case the minus-sign condition is satisfied if we choose the 
edge corresponding to the field X12 to come with a (—1). Then, denoting the black vertices 
1,2,3 from left to right and denoting the white vertices 1,2,3 from bottom to top, we can 
write down the Kasteleyn matrix: 

/I z w'^ \ 
K'^^'{z,w) =11-2 1 . (2.2) 
\w 1 z-^ J 
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Figure 2: The brane tiling (dimer model) for dPi. 



2.3 Characteristic polynomial and the toric cone 

Rather than Kij itself, the object we really need is the determinant of K, which is the 
characteristic polynomial P{z,w). Since P{z,w) is a polynomial in two variables, each 
monomial z°'w^ gives a point (a, b) in a two-dimensional lattice. These points make up a 
convex polygon called the Newton polygon, which in turn determines a toric variety. For 
details on the map between the toric variety and the toric cone we refer the reader to [8]. 
This completes the circle of correspondences. The beautiful fact is that this toric variety is 
precisely the singular Calabi-Yau cone to whose tip the stack of D3-branes is transverse by 
construction. 

We embed the Newton polygon into the z = 1 plan^ in Z'^, and denote the boundary 
points of the k-gon hj Vq, . . . , Vk-i, ordered anti-clockwise around the polygon. Set Vf, = Vq. 
Linear combinations of the vectors Vi form the toric cone C of the variety; for details see 
again for instance [8]. For our long- running example, the del Pezzo quiver, the characteristic 
polynomial is 

P'^^' {z, w) = det Kff' = -4 + z~^ + w"^ + z + zw, (2.3) 

and the toric cone is shown in Figure [3l 

In addition to the toric cone C, we also need to work with its dual cone C*. The edges of 
the dual toric cone are generated by the interior pointing normal vectors to the exterior faces 
of the toric cone, normalized so that their magnitude is the smallest that can be achieved 
with integer components. In other words, 

Wi^i+i ocViX Vi+i. (2.4) 

^ There is nothing special about the 2 = 1 plane: one can show that for a toric Calabi-Yau manifold 
the endpoints of the generators of the toric cone lie in a 2-plane, and without changing the variety one can 
perform a S'L(3,Z) transformation to map this plane to the z ~ I plane. Also note our abuse of notation: 
we label the coordinates on Z'^ by (x, y, z), but this z is not related to the dummy variable z in the Kasteleyn 
matrix. As we will not utilize the Kasteleyn matrix beyond this point, we hope this does not cause any 
confusion. 
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z = l 
plane 




t>l = (-1,0,1), t>2= (0,-1,1), t'3= (1,1,1) "12 = (1,1,1) "23 = (-2, 1,1) 

J)4 = (0,1,1), (0,0,1) "34= (0,-1,1) "41 = (1,-1,1) 

Figure 3: Left: The toric cone of dPi. Right: The dual cone. The filled-in circles are the 
edge generators normal to the faces of the toric cone, and the empty circles are the remaining 
generators. Note that for this example the generators of the toric cone lie on a plane; this 
does not hold for a general toric variety, but only for Calabi-Yau ones. 



To find the other generators of the dual toric cone, one simply needs to include the minimal 
set of interior points so that all the integral points in the dual toric cone are spanned by the 
generators. This is most easily illustrated by an example: in Figure [3] we plot the generators 
of C* for the del Pezzo quiver, with filled circles for the edge generators, and empty circles 
for the interior points needed to complete the set of generators of the dual cone. The crucial 
fact that enables the counting of operators from toric data is that every integer point in the 
dual toric cone corresponds to a unique single-trace BPS operator [1]. 

Returning to our example, in Figure H] we show the first three levels of points for the 
del Pezzo quiver. The single point on level is naturally the identity operator. On level 1 
the reader should compare the {x, y)-coordinates of the points with the winding numbers of 
operators in Table [TJ they match perfectly. Clearly the higher levels will match as well, due 
to the conic structure of the construction. 

Let us still fiesh out this correspondence between points in the cone and operators by 
writing down the operator corresponding to point (x, y, z). This operator is given simply by 

{x,y,z) ^ Ai^...Ai^, with Wi^ + ... + Wi^ = {x,y,z), (2.5) 

where the A's are the independent minimal loops defined in Section 12.11 and Wi is the gener- 
ator of the dual toric cone corresponding to the minimal loop Ai. Thus the correspondence 
above states that to construct the operator corresponding to the point {x,y,z), one has 
to choose a set of generators (minimal loops) such that the sum of these generators gives 
{x,y, z). There are generally many ways of choosing the loops so that this condition is sat- 
isfied, but the beauty of the construction is that all these operators will belong to the same 
equivalence class under the F-term constraints, which makes it possible to say that each 
point corresponds to a unique operator in the quiver. 
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z=2 



y 



z=0 ' • • • 



• • » w 

X • • • • 



Level Level I Level 2 



Figure 4: The first three levels of the dual toric cone of dPi, i.e., the slices z = 0,1, 2. Here 
the level counts the number of minimal loops. 



2.3.1 Z-minimization and the Reeb vector 

The final element we need is the Reeb vector b of the toric variety. Its significance is that 
it generates an isometry in the bulk geometry corresponding to the R-charge of the brane 
world-volume SCFT. Thus, in the dual toric cone, the Reeb vector is normal to planes 
containing operators of constant R. The directions generating an equal R-charge plane are 
then the remaining flavor charges [6], defined up to a residual SL{2, Z). As shown in [9], the 
Reeb vector b = {bi, b2, 3) minimizes the volume cut from C* by a plane 2(6, n) = 1. This 
is equivalent to minimizing the quantity Z: 

where {vi,Vj,Vk) is the 3x3 determinant constructed from the three column vectors. Z 
essentially corresponds to the Einstein-Hilbert action for a metric h on X, and thus mini- 
mizing it, i.e., finding the critical points, corresponds to finding the Sasaki-Einstein metrics 
for this action. In the above we set 63 = 3, which holds in the 5*17(3, Z) frame in which the 
Newton polygon is embedded in the plane z = 1. Thus finding the Reeb vector is reduced to 
an algebraic problem of minimizing a function of two variables. For dPi, the minimization 
yields the Reeb vector 

ft'^^i = (0,4- v^,3). (2.7) 

If we assemble the three U{1) charges of the BPS operator into a vector n G C*, then the 
value of the R-charge (and conformal dimension A) is given by its scalar product with the 
Reeb vector: 

A{n) = -R{n) = {b, n) . (2.8) 

Exemplary values of the R-charge for integer points in the dual toric cone of dPi are tabulated 
below. Notice that in general the R-charge of an operator need not be integer-quantized or 
even rational |10j . 

2.4 Toric cone as the phase space of dual giant gravitons 

So far we have concentrated on the description of the half-BPS spectrum of mesonic scalar 
chiral primary operators Timesonic in A/" = 1 SCFT quiver gauge theories. Our starting point 
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Level 


Vector 


Conformal Dimension 





(0,0,0) 





1 


(a,l,l) 


7- 


1 


(a,0,l) 


3 


1 


(a,-l,l) 


-1 + 


2 


(a,2,2) 


2(7- v^) 


2 


(a,l,2) 


10- 


2 


(a,0,2) 


6 


2 


(a,-l,2) 


2 + yT3 


2 


(a,-2,2) 


2(-l + v^) 



Table 2: Conformal dimensions of charge vectors n in the dual toric cone of dPi. The level 
denotes the number of minimal loops. The R-charge is independent of the first coordinate 
of n because bi = 0. 

was the mathematical description of toric cones Ai over Sasaki-Einstein manifolds X probed 
by D3-branes. However, the half-BPS sector of the gauge theory Hilbert space may also be 
analyzed in terms of (dual) giant gravitons. In this analysis, the machinery of toric geometry 
is as useful as it is on the gauge theory side. 

We recall that BPS mesonic sectors in SCFTs allow different physical interpretations 
depending on the amount of R-charge carried by the states. For states carrying order 
charge, these are interpreted as pointlike rotating gravitons; when the charge becomes of 
order N, they can also describe giant gravitons or dual giants. Classically, the first set 
corresponds to rotating D3-branes wrapping homologically trivial 3-cycles in the cone Ai 
over the base X; the second set describes rotating D3-branes wrapping a 3-sphere in AdSs. 
This last approach was adopted in [1], where the authors showed that Ai has a natural 
interpretation as the phase space of the dual giant. 

Using geometric quantization, the Hilbert space for the dual giant Tidg was found to be 
given by the space of holomorphic normalizable functions on A4. Since Ai is toric, it was 
further shown that this space is spanned by the elements of the dual toric cone C*. Thus, 
there exists an isomorphism between Tidg and the dual toric cone C*. 

Since dual giants are mutually BPS, the A^-dual giant states can be described by N 
indistinguishable quantum particle states. Hence the A^-dual giant Hilbert space is simply 
the A^-th symmetric tensor product 

n = Sym^7^,, . (2.9) 

This fact allows to compute the partition functions of these systems. We will be interested 
in the grand canonical partition function for multi-dual states, which was found in [1] and 
agrees with the grand canonical partition function of the CFT presented in [1]. We will be 
analyzing these partition functions in Section HI and interpreting the results both in terms of 
typical mesonic CFT operators, and in terms of typical configurations of dual giant gravitons. 
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Isomorphism of Hilbert spaces: Besides the isomorphism of Hilbert spaces described 
above, there exists a further one between the Hilbert space of giant gravitons and Ti [2]. 
This can be seen as follows. 

The classical moduli space of giant gravitons is given by the set of polynomials of degree 
N in A4 subject to a set of constraints [2] 

P{zi, Z2, Z3) = c + CiZi + CijZiZj H = ^ cizi . (2.10) 

Following [11], the quantum Hilbert space for such giants Tig is spanned by the states 

Ic/i, C/2, ...,C7^). (2.11) 

These states are holomorphic polynomials of degree N over the classical moduli space and 
are symmetric in the {c/.}. Thus, Tig is the symmetric product 

= Sym (|C7,) ® |c,,) ® ■ ■ ■ ® \ci^)) , (2.12) 

where each Ic/.) represents an holomorphic function over the cone M. with base X. This 
establishes the isomorphism between the Hilbert spaces of giant gravitons and dual giants. 

We want to emphasize that this last isomorphism is captured by the combinatorial iden- 
tity, a cornerstone of the plethystic program: 

gN{t,M)=gi{t,M''/SN), (2.13) 

where 5'Ar(t, M) stands for the partition function of the gauge theory of D3-branes probing 
M., whereas (71 (t, Ai'^/S^) stands for the partition function of a single D3-brane in the 
symmetric product M.^ ISn- 

The full discussion is summarized by the equality : 

Tig = Sym^TYdg = ^mcsonic • (2-14) 

We conclude that the partition function 5'Ar(t, M) counts half-BPS mesonic scalar chiral 
primary operators, and giant graviton states carrying a certain R-charge, and dual giant 
graviton states carrying the same R-charge. 



3 Counting multi-trace operators 

In the previous section we learned how each point n in the dual toric cone corresponds to 
an operator in the quiver theory via (12. 5p . These operators are not gauge invariant as each 
field transforms in a non-trivial representation under two of the gauge groups. In order to 
create gauge invariant operators we need to contract appropriately the indices on the fields 
Xi. This can be done by inserting traces at appropriate places in the operator El 

^ By constructing gauge invariant operators using traces we are picking a preferred basis that is convenient 
for us and which is related to multi-particle states in AdS space. Another possible choice for constructing 
such operators would be to use determinants and sub-determinants which are related to D-brancs in AdS 
space. 
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The simplest way of making the operator in fl2.5p gauge invariant is to take the trace of 
the entire quantity. This leads to a single-trace operator, which corresponds canonically to 
a point in the dual toric cone. As mentioned below eq. fl2.5p . this correspondence is bijective 
and extends to all the integer points in C*. Therefore every integer-coordinate point in the 
dual cone is uniquely associated with a single-trace gauge invariant operator. However, we 
can also create a gauge invariant operator by distributing several traces in the product (12.51) 
to produce multi-trace operators of the form (12. ip . This may clearly be done in a multitude 
of ways. In this section we will compute the number of such inequivalent trace structures. 
This is important because for large charges the corresponding field theory states are expected 
to be dual to black hole microstates, and our computation of the degeneracy of operators for 
any triplet of charges provides a prediction for the entropy of these black holes. We return 
to this in more detail in Section [5l 

Before proceeding, we should explain how our counting differs from the one done in the 
plethystic program [H [3]. There, the authors computed the entropy of mesonic operators 
of a given R-charge and provided an implicit expression-the refined plethystic exponential- 
counting operators at any given point in the dual cone (meaning operators of a given R-charge 
and flavor charges ni and 722). In this paper, we pursue a different strategy, which involves 
introducing a new slicing of the dual toric cone. This method affords an elegant geometric 
interpretation, and, crucially, allows us to find the asymptotics of the refined counting, which 
was not done in the plethystic program. As one expects the flavor charges of black holes to 
be measurable to classical observers, it is the refined counting carried out in this paper that 
is expected to match the entropies of the dual black holes. 

There is a simple way of stating our counting problem: since each trace component in 
the expression (12.11) is a single-trace operator, it corresponds to a point in the dual toric 
cone. Thus, computing the number of ways of distributing traces is equivalent to finding 
the number of partitions of the vector n into components: n = rhj. The asymptotics of 
this type of counting problems can be tracked using a theorem due to Meinardus [12]; our 
application of Meinardus' theorem will follow that of [1] . We shall proceed in the following 
steps: First, we count the number of partitions for all vectors on a fixed plane; this is 
equivalent to counting multi-trace operators with a single fixed charge, say R-charge. We 
then go on to count the vectors on this plane individually; this is the refinement, wherein 
a particular triplet of t/(l)-charges is fixed. We carry out the computation in this way 
because, as we shall see in the ensuing arguments, this method leads to a handy geometric 
construction involving centers of gravity of certain pyramids. 

3.1 Step 1: Partitioning co-planar vectors 

As the first step we shall reproduce the non-refined plethystic counting, and use the Meinar- 
dus theorem [12] to find the asymptotic number of multi-trace operators on any plane of 
rational slope containing the point n in the dual toric cone. Let us pick such a plane and 
denote it by p. Then consider an auxiliary family {pm} of planes parallel to p, such that (1) 
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each point in the dual cone is contained in one of the planes, and (2) the planes are evenly 
placed, i.e., the distance between any two neighboring planes is fixed. In this auxiliary family 
we index the planes from the tip of the cone, and take pn = p to be the plane that contains 
the point n. 

To illustrate this, we display a two-dimensional analogue in Figure [5l here the cone is 
chosen to be the positive quadrant of Z^, and parallel planes are drawn as parallel lines. 
Next we wish to find the sum of the number of different vector partitions (which we will 
call the multiplicity) of all the points on the plane p„. This is the first step on the way to 
finding the multiplicity of the point n. As is clear from the above setup and Figure [5], this 
problem is equivalent to enumerating the sets of vectors {mj} such that '^j'nij is on the 
plane Pn- Let rrij index the plane that contains rhj, counting from the tip of the cone. Then 
this condition is equivalent to 

rrij = n. (3.1) 

j 

Let us denote the number of integer coordinate points on the plane Pm by am- Then the 
problem is equivalent to the problem of partitioning an integer n in integer components rrij, 
with the added complication that at level m there are inequivalent summands. This is a 
standard partition problem and the asymptotics is studied by Meinardus [12]; his theorem 
states that asymptotically the number of such partitions is given by 

S{n) =log{# partitions of n} = (^AnT {a + 1) ({a + + O {log n) . (3.2) 

Here A is the residue at the rightmost pole, located at s = a, of the Dirichlet series D{s) 
given by 

oo 

D{s)^J2-s- (3-3) 

m=l 

Therefore, to find the multiplicities we need to know the coefficients a^, i-e., the number of 
(integer coordinate) points on plane Pm- In the plethystic program, the integer corre- 
sponds to the number of single-trace operators at level (distance) m, and the multiplicities 
we are after are integers dj such that 

oo oo 

JJ(l-r)-'*™ = J^rfjt^ (3.4) 

m=l j=0 

Here the exponents of t parameterize the distance of a plane parallel to p from the tip of 
the dual cone. The quantities dj are the numbers of multi-trace operators composed of the 
single trace ones. We will revisit this in Section HI 

Example: Let us see how this works for the two dimensional case portrayed in Figure [5l 
where we have chosen the lines (planes) to have slope —2 and the numbers am sue given 
by = 1, 1, 2, 2, 3, . . . —>■ 1 + [m/2]. We see that indeed, ao, ai, 02, . . . correspond to the 



12 



Figure 5: A two-dimensional partition example with one possible family of parallel lines with 
slope —2. The dark lattice points are the coordinate points and the numbers label the lines: 
Pn means the n-th line with slope —2 away from the origin. We have specifically chosen to 
draw the line pio as a solid line while other lines Pm are dotted and for piQ, we have portrayed 
one possible partition of vectors summing up to the lattice points on it. 

number of points on the lines po,pi,p2, ■ ■ ■ Clearly, for a general slop^ the coefficients 
are given by 

Om = 1 + 



m 
a 



- + 0{ll (3.5) 
a 



where [x\ denotes the integer part of x. 

Plugging ttm into (13. 3p . we can evaluate the Dirichlet series to give 



D(s) = -C(s-1) + ---, (3.6) 
a 

from which we find the location of the rightmost pole at s = 2, with residue 1/a, i.e., a = 2 
and A= 1/a. The omitted terms give additional poles which are further to the left. Thus 
the Meinardus theorem yields 

S{Pn) = ^(^r(3) C(3))' + O(logn) . (3.7) 

Note that the quantity An'^ = n^/a in the Meinardus theorem is equal to twice the area 
of the triangle cut off from the cone by the Pn line. This will be true also in the three 
dimensional case: the coefficient An" will be proportional to the volume of the pyramid cut 
off from the cone by the plane pn- This will be of critical importance later; it allows us to 
sidestep intractable analytic computations via an elegant geometric construction. 



We assume a > 1; the opposite case clearly works similarly. 
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3.1.1 Explicit enumeration of the planar partitions 

Let us now generalize from the two-dimensional example to a (i- dimensional toric variety. 
For us the interesting dimension is d = 3, but working in a general dimension comes at no 
additional cost. 

We start with the coefficients am'- the number of integral points of the dual toric cone 
contained in the plane Pm- For large m, the number of integral points can be approximated 
by the {d— l)-dimensional area of the plane pm- By dimensional arguments, this area scales 
as ttm = Am'^~^, where .A is a constant specified by the orientation of the plane in the cone. 
This allows us to evaluate the Dirichlet series (13.31) as 

h corrections = AC,{s + 1 — d) + corrections . (3.8) 

771=1 

This has the rightmost pole at s = rf, with residue A, yielding the parameters a = d and 
A = A. 

To evaluate A, consider the number of integral points inside the pyramid cut off from 
the dual toric cone by the plane Pm- This can be approximated by the volume Vm of the 
pyramid. Again by dimensional arguments, for a d-dimensional cone this volume scales as 
Vm = Am"' for some A. Clearly the number of points in the plane pm can then be written as 

am = Vm- Vm-1 = A(m"' - (m - l)'^) ^ Arfrn'^-^ + ^(m'^-^) ^ (3.9) 

m 

from which we find the residue 

A^^^^. (3.10) 

Using these results we can plug into (13.21) and write the number of partitions of vectors on 
plane p as 

S{p) = ln(# partitions) = (vdT{d + 1) ({d + 1)) ^ + 0(log V) , (3.11) 

where V is the volume of the pyramid cut off from the cone by the plane p. We have dropped 
the index n in p„ because it is a function only of the particular plane in question, and not 
the additional family of auxiliary planes used in the construction. 

We will need an explicit equation for the volume V later in Section H] and will give the 
somewhat complicated expression there. For now we do not need such an explicit formula; 
it would even be a hindrance, as the analytic computations would quickly become hopelessly 
complicated, while simple geometric arguments can take us remarkably far. 

3.2 Step 2: Partitions of individual vectors 

Having found the partitions of all vectors lying on a given plane p, we now wish to use the 
results above to find the possible partitions of an individual vector n. We shall accomplish 
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this by considering an upper bound on the number of such partitions, and then showing that 
asymptotically this bound is saturated. 

Consider {pfi,s}, the set of all planes containing the point n. Here the subscript s indicates 
that the plane Pft^g is defined as being normal to s, as well as containing n. Note that for 
Ph,s to cut off a finite pyramid from the cone, we must have (s, Wa) > for all the edge 
generators Wa of the dual coneB One can see that this condition is equivalent to s living 
inside the toric conejf] Continuing the above notation, we denote the volume of the pyramid 
cut out by the plane Pfij as Vfij. 

Since the point n is merely a point in the plane Pfi^g, the multiplicity of n is bounded by 
the multiplicity of pn,s- This gives: 

^sin) < g5(pH..) = exp (^^^ {dVH,sT{d + l)C{d+l))^^^ . (3.12) 

Since this is true for any plane Pft^g, we can take a minimum over all the possible s on the 
right hand side of (13.121) . This minimum is well defined because s lives inside the toric cone, 
and on the faces of the cone Vfi^g diverges]^ Thus a minimum will exist inside the cone. 
Defining the function f{n) as 

/(n) = dmmVH,g, (3.13) 

we can write f l3.12p as 

eS(n) < (/(^)r(d + i)((d + i))dTT^ . (3.14) 

Now we wish to argue that in an asymptotic limit this bound becomes saturated, and there- 
fore (13.141) gives us the multiplicity we are after. 

It is straightforward to see which s minimizes the expression (I3.13P : by construction the 
volume minimizing plane, PH,smin^ to be tangent to surfaces of constant /, and therefor^ 

Smin = V/(f2), SO that f {H) = dV^^^j^^y (3.15) 

Thus, to any point n we can associate the volume minimizing plane 

p{n) = pH,g^^^ . (3.16) 

In Sec. 13.3.11 below we show that n is the center of gravity of the polygon p{n) fl C*. The 
uniqueness of the center of gravity then establishes that p{n) is injective. But this implies 



Since s and —s define the same plane, {s,Wa) < would work just as well. We choose > without loss 
of generality. 

^ Note: not the dual cone C*, but the toric cone C. 

^ This is because for s € dC, the constraining planes become parallel to faces of C*. Thus, they fail to 
contain finite volumes within C* so Vn.g diverges. 

^ To be precise, this argument only implies that Smin is proportional to the gradient V/, but since the 
normalization of s is immaterial, we can use this equation to fix the normalization. 
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that the point n maximizes the multiphcities on the plane p{n). This is because, by the 
injectivity of the map p, all other points n' on the plane must have different minimizing 
planes p{n') for which n' will be the center of gravity. Then, necessarily, \^(^') < since 
p{n') is volume- minimizing plane for n'. Thus, 

g5(n') < g5(p(n')) < g5(p(n)) _ f (/(^)r(rf + l)({d + 1))^^ • (3.17) 



This then implies that 



where A(p(n)) is the area of the volume minimizing plane passing through n. The left hand 
inequality simply uses the fact that S{n) > S{n') on the plane p{n) while e'^*-^^"^'' = e"^*^"'-*. 
The right hand inequality is fl3.14p . Now observe that the A{p{n)) scales as |f2|'^~^ where we 
are taking the Euclidean norm. Thus, taking logarithms of both sides of f lS.lSp we get 

S{p{n)) - C(log(|n|)) < S{n) < S{p{n)) (3.19) 

Now from (13.111) S{p{n)) scales as V~ ~ |n|*+T. Thus we see that when \n\ is large, 

gS(n) ^ g5(p(n)) ^ (/(^)r(rf + l)((d + 1))^^ , or 

S{n) ^ ^^(/(ri)r(d+l)C(rf+l))^ . (3.20) 

This is the main result of this section. We have only kept leading terms in all expressions; 
to accurately evaluate multiplicities near the boundaries of C* subleading terms will have 
to be kept too. Eq. (13.201) can be used to compute the number of multi-trace operators 
corresponding to any given triplet of charges, and thus gives a prediction for the entropy of 
the dual black holes on the gravity side. 

Z-minimization revisited: Recall from Sec. 12.3.11 that the Reeb vector b = (61,62? 3) is 
computed by minimizing Z, which is proportional to the volume cut from C* with a plane 
of the form 2(6, n) = 1 |4]. Notice that all such planes pass through the point n = (0, 0, |). 
The Reeb vector is therefore precisely the normal to the plane p((0, 0, |)) and indeed to all 
planes p(0, 0, z): 

6cx V/((0,0,;2)), (3.21) 

where the proportionality is set by requiring that 63 = 3. The argument above eq. (13.171) 
then establishes that every point on the 2-axis is the maximal entropy point on its own equal 
R-charge surface. Because these surfaces take the form 

(6, n) = b,x + 621/ + 3z = — , (3.22) 
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Figure 6: Equal entropy surfaces in the dual toric cone of C'^. They take the form of 
generalized hyperboloids (in the interior of the cone; see left) and hyperbolae (on the faces 
of the cone; see right.) The interior surfaces represent parametrically higher entropies. 



the point of maximal entropy of an equal R-charge surface is always given by 

(0,0,:^). (3.23) 

In the next subsection we generalize this to arbitrary constraint surfaces. The result (13.231) 
holds in the standard presentation of the dual toric cone, that is when the Newton polygon 
lives in the plane z = 1. 

Surfaces of constant entropy: From f l3.20p we see that the surfaces of constant entropy 
in the dual toric cone are given by level sets of the function f{n). From the definition of 
f{n), these surfaces are characterized by the property that at each point n on the surface, 
the plane tangent to the surface at n cuts out an equal volume from the cone. Surfaces with 
this property are generalized hyperboloids in the cone. In Figure [6] we plot these surfaces of 
constant entropy. 

It is also interesting to consider entropy on the faces of the cone. Although the preceding 
construction appears to assign zero entropy to these faces, this should taken as a statement 
that the entropy on the faces is parametrically smaller than inside the cone. A good example 
here is C^, where the interior of the cone corresponds to the |-BPS sector, i.e., operators 
composed of three fields {X,Y,Z}, while on the faces only two fields (say {X,Y}) are available. 
Evidently the entropy on the faces is parametrically lower but non-zero. The equal entropy 
curves on the faces take the form of hyperbolae. (See Figure El) 

3.3 Maximizing entropy 

In Section [32] we saw that the degeneracy of multi-trace operators with a specified rf-tuple of 
charges n is entirely determined by the behavior of one function, /(n), which was essentially 
given by the minimal volume that can be cut out from the dual toric cone with a. {d— l)-plane 
containing n. The normal to this minimizing plane is then given by the gradient V/(n), 
and it is very useful to think of V/(n) as the direction of the most rapid growth of entropy, 
starting from the point n. 
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A problem of interest is to maximize entropy subject to a specified constraint, such as 
a fixed R-charge R. In the following we consider constraints which are linear in the charge 
vector n. Such constraints trace [d — l)-hyperp lanes in (i-dimensions and take the form 

c,-(n) = -n, = 0, (3.24) 

where s is the normal of the constraint plane, and is proportional to the corresponding 
charge, with the proportionality constraint set by the normalization of s. For the special 
case of fixed R charge, the normal to the constraint planes is given by the Reeb vector h of 
the toric variety, and the constraint equation reads (fe, n) — |-R = 0. 

Maximizing entropy given a constraint cg{n) = is then accomplished using the standard 
Lagrange multiplier technique: the maximum is given by the point where the constraint plane 
is tangent to a surface of constant f{n) which is equivalent to 

V/(n) oc Vcs{n) = s. (3.25) 

Unfortunately, due to the complicated analytic expression for the volume of a pyramid 
cut from the cone, equation fl3.25|] is difficult to work with analytically. Luckily, a simple 
geometric observation reduces this problem to a very tractable one, familiar from freshman 
mechanics. 



3.3.1 Maximal entropy points as centers of gravity 

We shall now show that the point of maximal multiplicity on a given plane is given by the 
center of gravity of that plane. This provides a simple prescription for finding the maximal 
entropy given an arbitrary constraint. Using the AdS/CFT correspondence, this will give 
the maximal entropy of a black hole in AdS space with charges satisfying the constraint. In 
particular, this will allow us to find the entropy and fiavor charges of the most entropic black 
hole with R charge R. 

Consider a linear constraint cg{n) = 0, and denote the constraint plane by pg. Let the 
center of gravity of pgCi C* be Hq. Now, pick an arbitrary lin^ / contained in the plane 
pg that passes through the center of gravity no- This is portrayed in Figure [71 To find the 
center of gravity, we define coordinates on the plane pg. We do so by choosing one of the 
coordinates, x, to be along I, and the other, y, to be perpendicular to /. Then the center of 
gravity is defined by the condition 

dxdyy = 0, (3.26) 
where the integration is over the part of the plane pg contained in the dual toric cone, i.e., 

A = pgnc*. 

^ For clarity, we present this argument for the three-dimensional cone. However, the same logic will go 
go through in an arbitrary dimension d. 
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Figure 7: Left: A constraint plane slicing the dual toric cone. Right: A line / (in blue) 
through the center of gravity of the slice. The contour in red illustrates the change in the 
volume cut by the plane after it is tilted about I by the angle a. The toric variety used in 
the present example is the cone over dPi. 

Now, consider tilting the plane pg around the axis / by a small angle a, as illustrated in 
Figure (^^. Denoting this tilted plane by p'g, the change in volume cut from the cone can be 
written as 



where A' is the intersection of the tilted plane p'g and the dual cone C*. The approximation 
in the equation above requires a little bit of explanation, as one needs to take into account 
the change of region of integration from A' to A. However, it is easy to see that this change 
will contribute terms that are proportional to a, and combined with the sin a factor these 
corrections only show up at order (9(a^). The condition that pg is the volume minimizing 
plane, AV = 0, is therefore equivalent to setting the torque about the axis / to vanish as in 



However, since the line / is an arbitrary line passing through the center of gravity no, 
the above argument is valid for any such line, and therefore for any tilt about no- It follows 
that the constraint plane pg is the volume minimizing plane for the center of gravity no , and 
therefore, the arguments from Section 13.21 tell us that no must be the most entropic point 
on the plane pg. 

As a side note, the preceding argument affords a new characterization of the Reeb vector. 
In particular, b is the normal to the unique family of parallel planes Pn such that the center 
of gravity of pn fl C* is contained in the 2;-axis (x^-axis in d dimensions) This observation 
is applicable in the standard presentation of the toric fan of Ai, where the Newton polygon 
is contained in the plane z = 1. 

In some practical situations, finding the center of gravity of an infinitesimally thin slice 
(s, n) = c of the dual toric cone may still be harder than solving another auxiliary problem. 
Consider a family of parallel slices (constraints) of the form (s, n) = q, indexed by i. Because 
C* is a cone, all such slices are similar polygons, differing only by an overall scale. This 

^ With a slight abuse of notation, we use z to denote a coordinate both on the toric cone and the dual 
toric cone. 




(3.27) 
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Figure 8: A family of parallel planar constraints intersecting the dual toric cone. All con- 
straint slices are similar polygons, so their centers of gravity fall on a ray extending from 
the tip of the dual toric cone. Here the toric variety is dPi. The constraint slices are normal 
to the Reeb vector and represent equal R surfaces. The diagram was drawn in the standard 
5*^(3, Z) frame, so that the centers of gravity of equal R-charge slices fall directly above the 
tip of the cone. 

means that their centers of gravity all fall on the same straight line k{s), which is uniquely 
determined by the normal vector s of the constraint plane. The intersection of this line with 
the constraint (s, n) = c then yields the maximal degeneracy point on the plane. This is 
illustrated in Figure [H 

Now consider the pyramid Ag^c cut out from the dual toric cone by the plane of constraint 
{s,n) = c. Its center of gravity also evidently lives on the line k{s). We may now reverse 
the argument and find the line of maximal entropy, for this family of constraints given by 
(s, n) = Ci, by finding the center of gravity of the pyramid c. In particular, k{s) is the ray 
extending from the tip of the cone in the direction of the center of gravity of A^^ c- In this way 
the point of maximum entropy on the constraint (s, n) = c may be mechanistically recovered 
by hanging the pyramid by its tip using a piece of string. When the solid is balanced, the 
maximal entropy point will be directly below the tip. 

The latter observation leads to a highly peculiar, mechanistic characterization of the Reeb 
vector. Consider the dual toric cone in standard presentation and align the 2;-axis (x^-axis 
in d dimensions) with the vertical direction. Of all families of pyramids Ag^c, the family Ag ^ 
is the unique one whose centers of gravity will lie on the vertical axis, i.e., directly above 
the tip. In other words, if we chop the dual cone C* with a plane normal to b, the resulting 
pyramid will balance on its tip! 
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Figure 9: Two sets of constraints on the dual toric cone of C^. Level sets of the R-charge 
(normal to b = (1, 1)) and of the quantity x + 2y = |i? + y are given in blue and red, 
respectively. The midpoints of each set of line segments (bold dots) trace rays of fastest 
entropy growth for each set of constraints (dotted vectors). 

3.3.2 Examples 

We begin with the simplest possible example: and take the dual toric cone to be the first 
quadrant of Z^, which is not the standard presentation. Here the Reeb vector is given by 
b = (1,1), and thus the lines corresponding to the constraint of fixed R charge are normal 
to (1,1), i.e., they have slope —1. Here every equal R-charge constraint is a line segment 
normal to the Reeb vector b = (1,1). The center of gravity of each line segment is simply 
its midpoint, and the union of all midpoints forms the fastest growth ray k{b) = (1,1). 
This is drawn in blue in Fig. [91 The numeric equality between 6 = (1, 1) and k{b) is purely 
accidental^ and is not respected by SL{2, Z) transformations. One may wish to impose a 
different constraint, for example of the form x + 2y = const. The ray of maximal entropy is 
again given by the union of the midpoints of the line segments x + 2y = c contained in the 
dual toric cone, and we can read off k{{l, 2)) = (2, 1). This is illustrated in red in Figure [H 
As a next example, consider the simplest three-dimensional toric variety. In its simplest 
presentation, where both the toric and the dual toric cone are simply the first octant of the 
lattice Z^, the symmetry relating the three generators is evident and trivially k{b) = (1, 1, 1). 
For a less trivial example, consider the presentation of the toric data of in the conventions 
of Section [2l where we choose the SL{3, Z) frame so that the generators of the toric cone 
reside in the z = 1 plane. It is straightforward to compute the generators of the dual toric 
cone and the Reeb vector using the results from Section [2], and one gets 

u;i = (0,1,0), ^2 = (-1,-1,1), ^3 = (1,0,0), 6= (1,1,3). (3.28) 

Choosing the slicing corresponding to fixed R charge, i.e., the one corresponding to planes 
normal to the Reeb vector, we find that the center of gravity of the pyramid cut out from 

These two objects should not be compared, as 6 = (1, 1) hvcs in the toric cone while k{b) lives in the 
dual toric cone. 
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C* by the plane of R-charge R, x + y + 3z = is the point (0, 0, ^). Therefore the hne of 
maximal entropy at fixed R-charge is given by 

fc(6) = fc((l, 1,3)) = (0,0,1). (3.29) 

By the argument below eq. f l3.2ip . this holds for every toric variety in the standard pre- 
sentation. We remind the reader that k{s) is to be understood as specifying a ray, so its 
normalization is immaterial. The ray intersects the equal R surface at 

{Q,0,^) = ^{w^+W2 + Ws), (3.30) 

in agreement with eq. fl3.23p . The form of (13.301) . symmetric in the generators Wi, W2, w^, is 
a direct consequence of the fact that the three adjoint scalars X, Y, Z have equal R-charge. 
As exemplified here, all algorithms presented in the current section respect the SL{3, Z) 
symmetry. 

Finally, we treat the case of dPi, which was used for illustration throughout Section O 
Again consider the family of equal R-charge constraints, which in the (standard) SL{3, Z) 
frame of Figure [3] take the form 

(4 - Vl3)y + 3z = ^Ro . (3.31) 

The center of gravity of the pyramid cut out from C* with the plane (13.311) again resides at 
(0,0, ^^); indeed, this will hold for all pyramids formed by equal R-charge surfaces in the 
standard presentation. The most entropic charge vector for a fixed R-charge Rq takes the 
form (13:231) . 

The appearance of an irrational number in the Reeb vector of dPi, eq. (12. 7p . implies that 
an exact specification of the R-charge traces a line segment and not a plane intersecting C*. 
Namely, all the integer points in the dual toric cone that line within a plane perpendicular 
to the Reeb vector are actually coUinear. In detail, writing down 

Ro = Ri + R2V13, (3.32) 

the locus of equal R-charge takes the form 

Ay + 3z = 

-y = ^R2, (3.33) 

which is a line in three-space. However, such considerations would only be of practical import 
if one were capable of making exact R-charge (conformal dimension) measurements, so as to 
extract from one measurement two integers -Ri, i?2- At any finite precision e, the relevant 
constraint to impose instead is 

Ro-e<R<Ro + e. (3.34) 
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In this case, the most entropic point is given by the center of mass of the thin shell fl3.34p and 
is given by (13.23p . For completeness, we present in Appendix |A] an algorithm for calculating 
the multiplicity of operators matching a specified irrational R-charge exactly (analogously 
to eq. fl3.1ip ). This generalizes the construction of Sec. 13.11 

3.4 Dual giant interpretation 

Recalling the construction reviewed in Sec. 12. 4[ wherein C* is seen as the phase space 
of a dual giant graviton propagating in the AdSs x X background, one can also interpret 
the present results in terms of dual giant configurations. Dual giant gravitons are point- 
like objects from the viewpoint of the Sasaki-Einstein manifold X and carry momentum 
P = (P(^^, P^g) along the directions of the maximal torus of X. This momentum is 
quantized and takes integer values falling in the dual toric cone C*. In this way, the conserved 
momentum of each dual giant graviton P corresponds to a charge vector m comprising the 
three U{1) (flavor) charges studied above. Likewise, we may think of a vector s* G C as 
specifying a particular cycle ag in the maximal torus, and a planar surface 

{s, P) = Pg (3.35) 

consists of all states whose component of momentum along ag is fixed at Pg. The result of 
Sec. 13.11 eq. (13. lip , gives the total number of such configurations. 

Analogously, Sec. 13.21 gives a way of enumerating the configurations of many dual giants 
which are characterized by a given triple of total momenta P = {P^p^ , P<^2, P^g). Following the 
procedure, eq. (13.150 assigns to P a special cycle corresponding to V/(P) G C. Eq. (I3.20p 
then shows that fixing the total momentum along that cycle becomes in the asymptotic 
regime equivalent to fixing the momentum vector at P. The reason for that is reminiscent 
of the equipartition theorem and may be traced directly to the central limit theorem. In 
particular, the orientation of the cycle crv/(P) maximal torus of X is such that mo- 

mentum running along it splits among the fundamental cycles of the torus precisely in the 
ratio given by P. 

4 Universal distribution of trace factors 

In the previous section we found the multiplicity of operators for any given triple of charges 
n, comprising (linear combinations of) the R-charge and the two flavor charges. Further, 
we gave an algorithm which selects the combination of flavor charges n that maximizes this 
multiplicity given any linear constraint on the charges. In this section we will show that 
when \n\ is large, most of these operators share a universal trace structure. 

In [13], the structure of multi-trace operators in the half-BPS sector of A/" = 4 SYM 
theory was treated using canonical partition function techniques. The problem of counting 
BPS multi-trace operators in an A/" = 1 SCFT living on a finite stack of D-branes transverse 
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to the tip of a Calabi-Yau cone was likewise treated in [T] using the plethystic program. In 
the following we shall proceed along similar lines in order to extract statistical information 
about the distribution of traces in a multi-trace operator with charge vector n. 

As explained in Section [3l every multi-trace operator of charge vector n is specified by a 
vector partition 

^mj=n, (4.1) 

3 

where each summand rhj e C* corresponds to one trace factor, unique up to F-term relations. 
The exact correspondence between flavor charges, winding paths in the dimer model, and 
fields entering the given trace factor, was described in Section [21 Introducing a triple of 
generalized temperature^^ /3 = (A, /32, /^s) and a chemical potential a , we may write down 
a grand canonical partition function for the system [H H] 

Z = J] (1 - e"-'^^'"''>y' = n • (4.2) 

This partition function has been exhibited in [1] as counting dual giant graviton states moving 
on a toric variety Ai, and concurrently in [1] as a generating function for mesonic multi- 
trace operators in the associated field theory. In the following, we first treat the case where 
a = so that every summand m (trace factor of type m) is counted with the same weight. 
Indeed, each factor = (l — e"^^'™"-*) allows for an arbitrary number of repetitions of 
the summand m in a partition (multi-trace operator). Thus, in setting a = 0, we are going 
to let the ensemble choose its optimal number of trace factors M. The case where M is fixed 
by hand with a non-trivial chemical potential is treated later in the section. 

4.1 Ensemble with unrestricted number of traces 

With a = 0, the partition function (14. 2 p blows up at G C*. This corresponds to the onset of 
Bose- Einstein condensation at the trivial point. However, the tip of the cone corresponds to 
a trivial factor in the structure of a multi-trace operator, and as such should not be counted 
in the total number of factors in a multi-trace operator. Therefore, the correct partition 
function to be analyzed is: 

Z= H n ^rn. (4.3) 

mGC*\0 meC*\0 

For a given charge vector n, the Lagrange multipliers (3 are set by requiring: 

i^i) - 

These temperatures are not physical 



They are Lagrange multipHers introduced for ease of calculation. 
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Figure 10: Illustration of the natural grading associated to the point (ni, n2, n^) for C^. The 
planes with dashed lines are normal to /5, and for comparison we have included the level 1 
surface containing Tr X, Tr Y and Tr Z with dotted lines. 



Because eq. (14. 3p factorizes, we can read off the mean and standard deviation in the popu- 
lation of the trace factor m: 

(n^,) = -^^l^ = (e(^'"^)-ir^ (4.5) 

The right hand side is manifestly independent of z. Indeed, fl4.4|] counts average charge while 
(14. 5 p counts the average number of traces; their relation is 

(rii) = ^ mi{n,7^) . (4.6) 

meC*\0 

The standard deviation in Urn reads: 



(4.7) 

so that the standard deviation always exceeds the mean, most rampantly away from the tip 
of the cone: 

a{n,^)/{nrn) = = (l + {n^r')'^' > 1 . (4.8) 

This is the same relation as that found in [13] . We defer a full discussion of the structure of 
heavy multi-trace operators until the next subsection. 
Eq. (14. 5p shows that the planar surfaces 

0,m) = y (4.9) 

share the same occupation numbers. This induces a natural grading of the dual cone C* by 
parallel planes with the normal P; we illustrate this in Figure [TOl It is easy to see that these 
planes are the same as the planes pm utilized in Sections 13.11 13.21 to evaluate the number 
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of vector partitions of a given vector n@ Consequently /3 || V/(ri) and we may recast the 
populations at sites m in terms of just one temperature 13: 

(n^) = -7^4—-^^ = (e^(-^/("^)) - 1)- , (4.10) 

where we have set 

/3 = /5V/(n). (4.11) 
With this form of (3 substituted in eq. (14.31) . the condition setting (3 reads: 

V/(n)) = -^^. (4.12) 

For large |n|, we can use a continuum approximation. A thin shell contained between 
two planes normal to V/(ri) given by 

t < {m,Vf{n)) <t + dt (4.13) 

is seen from eq. (13. 9p to contain 

a^dt = ^dt (4.14) 

lattice points. In order to write an explicit expression for Vf in d dimensions, we shall need 
further notation. Write the dual toric cone C* as a union of F simplicial cones a* and denote 
the d generators of a* by uJj, j = 1, . . . , d. Each of these is a d-vector, so {w^)jk is a. d x d 



matrix. In terms of the notation of Sec. 12. 3^ which is suitable to ci = 3, one choice of the F 
simplicial cones leads to: 

Wl = Wk-1,0 

iJj^ = Wi_i^i (4.15) 

Here i ranges from to F = k — 2, where the Newton polygon is a k-gon. The volume of 
the dual toric cone bound by the plane {s, x) = t with normal s is then given by 

^. = ^E4^^^. (4.16) 
It follows that the thin shell f l4.13p contains 

F 



det 



t'^~^dt = B{n)t'^~^dt (4.17) 



(^-i)'trn-=i(^;'W(ri)) 



This follows from the observation that the counting procedure in Sections 13.11 13.21 did not distinguish 
among charge vectors lying on the same plane p„j. Thus, the planes p,„ comprise vectors of equal multiplicity, 
as do the planes (|4.9p . so the two families must coincide. 
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m = {n^vmYB{n) = ^";;:\7 ri'..... • (4.19) 

i=i iii=i( 



lattice points. Recalling from the definition (13.151) that f[n) = Vtd ioi t = (n, V/(n)), we 
see that B{n) may be related to f{n) as 

B{n) = ^ = -JB— , (4.18) 
{n, Vf{n)r 

that is, 

( n, Wfin)Y y. \detWi 

Now, each of the points in the thin shell has mean population (exp/9t — 1)~^, and 
contributes t to the right hand side of eq. (14.121) . as a continuous version of (14.61) . Therefore, 
we may approximate the condition (14.121) for /3 by an integral: 

-^^dt={n,Vf{n)), (4.20) 

where the negligible contribution for large t was added to the integral to facilitate explicit 
evaluation. Evaluating the integral, we get 

^ Bin)Tid + l)C{d+l) 
(n, V/(n)) 

which, using eqs. (13.201) and (I4.19p . yields an identity: 

d+1 (S, V/(n)) 

Eq. (14.211) is in agreement with the well known case d = 1. There f{n) = n so that B{n) = 1 
and = ({2)/n, as in eq. (3.7) of [11]. In the case of primary interest, d = 3, in the 
notation of Sec. 12.31 /3 takes the form: 

V30(n, V/(n))(V/(n),«Jfc_i,o) (V/(n), ^._i,) (V/(n), ^.,+i) 
4.1.1 Structure of heavy multi-trace operators 

Eqs. (14.101 14.171 14.211) allow us to uncover a wealth of information about the structure of a 
typical large charge multi-trace operator. We list such observations below. 
The average number of trace factors in the operator is given by: 

, .^x r t'^^^ , Bin)Tid)C{d) 
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Now, using the relation between f{n) and B[n) from fl4.19p . we therefore have 

B(n)T(d)C(d) , . ^ ,,, d 

(M(n))o = ^ ^ ^(n,V/(n))^+i 

(5(n)r(rf+l)C(c/+l))^ 



d+l 



In the first hne we substituted for (3~'^ from eq. (14.211) and in the second for (n, V/(n)) from 
eq. (14.191) . The average number of traces is directly proportional to the entropy, with the 
constant of proportionality determined solely by the dimension of the toric variety d. 

The mean multiplicity (n^) of the m-trace in the large multi-trace operator, eq. (I4.10p . 
quickly drops to as |m| increases. The only traces whose occupations are of the order of 1 
or greater are those satisfying 

/3(m, V/(ri)) < 0{1). (4.26) 

By eq. (14. 8p . the same relation with a <S sign defines the regime where the mean accurately 
predicts a non-zero occupancy. Substituting for P from eqs. (14.211 1^^221) . we have: 



m, V/(n)) < _ v^..^^^^^ (^, V/(n)) 



(n,V/(n)) - dSin) ^ ' ^^^^^^ ~ lr(rf +\) C(rf + 1) i?(n) " ^^'^^^ 

This scaling behavior motivates the introduction of an auxiliary variable 

y{m)=(3{m,Vf{n)). (4.28) 

In comparison with the variable t of eq. (14.131) . y is rescaled hj y = (3t. It is clear that mean 
occupation numbers and all other quantities which are extractable from Z{j3) are functions 
of y alone. It is convenient to introduce the quantity 



<y) = / • (4-29) 




By the same logic which led to eqs. (14.201 14.240 . the quantity x{y) gives that fraction of 
trace factors making up a heavy multi-trace operator (given by (I4.24p ) which have y < y. 
The graph of x{y) is therefore a convenient presentation of all the universal features of a 
heavy multi-trace mesonic BPS operator. Importantly, x{y) depends only on the dimension 
of the toric variety d, which motivates the notation Xd{y)- This leads to the interesting 
conclusion that in all world-volume SCFTs on D3-branes transverse to tips of toric Calabi- 
Yau threefolds, all heavy BPS operators are described by the same curve xs{y). The curves 
X2{y), Xs{y), and X/i{y) are plotted in Figure [TTJ 
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Figure 11: Left to right: X2{y), Xs^y), x^i^y). 



The graphs Xd{y) are interpreted as follows (we leave a speculative interpretation to 
Appendix [C]). Suppose one is interested in the structure of a heavy multi-trace mesonic BPS 
operator with charge vector n. One begins by computing B{n) from eq. 04.171) and /3 from 
eq. (14.211) . The number of trace factors in equivalence classes fn satisfying 



Here we have differentiated eq. (14.291) and used eq. (I4.24p with y = [3t. Using eq. (14.171) one 
may also evaluate the average number of trace factors corresponding to an individual site 
m. This is given by 



Curiously, mean occupations per site depend only on the orientation of the vector n in C* . 

Since the quantity Xd{y) stands for the fraction of trace factors with y < y, it varies 
between and 1. In particular, 1 — x{y) is the fraction of operators with y exceeding y. 
From the fact that x{y) approaches 1 very rapidly for y of the order of unity one sees that 
almost all summands satisfy |m| oc The characteristic scale on which the occupations 

(I4.32P vary is thereby set by the quantity |n|^. 

Recalling eq. (I4.26p . the sites with occupations of the order of 1 or greater correspond to 
the section of the graph with y <1. These are also the sites for which the mean occupancies 
(I4.10[ I4.32P are good predictors for actual occupation numbers. Since ^3(1) ~ 0.15, we see 
that approximately 85% of the trace factors entering a heavy operator are randomly selected 
from a large volume in C* . The size of that volume is such that y > 1 but still of the order 
of unity. In particular, X3(10) ~ 0.998, which means that on average only 1 in 500 trace 
factors has (m, V/(n)) > 10 (n, V/(n))d+T. 

Denote the inverse of Xdiy) as yd{x), so that yd{xd{y)) = y- Consider the differential 
segment {x, x + dx). It stands for {M{n))odx traces, each contributing (3~^yd{x) to the 
right hand side of eq. (I4.20p . Therefore, yd{,x) dx is proportional to that fraction of the total 



f5-^y < (m, V/(n)) < [3-\y + dy) 



(4.30) 



is given by 




(4.31) 




(4.32) 
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charge of the heavy operator, which is carried by the traces with yd{x) < y < yd{x + dx). In 
particular, the charge carried by the trace operators with y not exceeding y is given by 

r^div) ry yd~i 
R'{y <y) oc yd{x) dx oc y— dy oc Xd+i{y) ■ (4.33) 

Here we performed a simple change of variable using the definition of Xd{y), eq. (14.29^ . 
Eq. (I4.33P holds for any charge of interest R' because for any linear combination of the d 
U{1) charges the total charge grows proportionately as one includes factors of increasing y. In 
conclusion, the function Xd+i (y) gives the fraction of the total charge of the operator residing 
in trace factors with y < y. As an example, ^4(1) ~ 3%, which means that only 3% of the 
total charge of a heavy operator is carried by trace factors whose expected occupancies are of 
the order unity or greater. The remainder is contained in trace factors which are randomly 
pooled from the volume (n, V/(ri))*+i < /3~^y < n (m, V/(ri))*+i , where k > 1. 



4.2 Ensemble with a fixed number of traces 

We now turn to the case where a < 0. The condition (14. 4p is supplemented by a condition 
fixing the number of traces M. In keeping with the standard language of the grand canonical 
ensemble, we shall alternately refer to trace factors as the particles of the system. The 
condition sets the parameter a: 

M = ^ . (4.34) 

oa 

Here Z is again as given by equation (14. 2 p and includes G C*. The results (14. 5114. 8 p are 
altered simply by substituting: 

0, m) 0, m) - a (4.35) 

and as a consequence, the planar surfaces (14.90 continue to comprise points of equal occu- 
pations. One consequence of that is that the validity of Sec. 13.3.11 will also extend to the 
present section: on each plane p intersecting C*, the charge vector of maximal entropy in 
the M-trace sector is the center of gravity of {p fl C*}. Furthermore, the argument which 
motivated (14. lip also continues to apply. Thus, we shall keep on referring to the scalar 

— * — » 

quantity /3 under the identification P = pVf^n). 

Proceeding as in Sec. 14.11 we approximate eqs. (14.41 14.34p by integrals: 

where gv{e°') stand for Bose-Einstein, or polylogarithm functions. We recall the definition 



00 ^ 

to be gy{z) := 1^- 

k=l 
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Using eq. (14.191) to eliminate (3 from the both equations, we get: 

^ ^ B^n)mU^') V/(,1))^ ^ ^ ^ (4,38) 

(B(,1) r(d + 1) gj+ile")) (r(d + 1) gj+ile")) 

Eq. fl4.38p allows one to set a independently of [3. However, the coefficient of f{n)'^ in 
fl4.38p is bounded from above by 

(r(<i+i)c(rf +!))"« 

since < < 1. For the case of most interest, c? = 3, we have 

fi:(3)^ 0.59. (4.40) 

This quantity marks the critical number of trace factors for the onset of Bose-Einstein 
condensation. Accordingly, we define 

M,Ht(n) = «(d)/(n)^ (4.41) 

as the critical number of traces for a given charge vector n. Note that in the a = case, 
wherein we let the ensemble choose its number of traces, we had 

(M(n))o = M„it(^), (4.42) 
as can be easily verified from eq. (14.251) . For M ^ Merit-, we distinguish the following cases: 



4.2.1 The classical phase 

For M < Merit (^), ('T'o) is a finite number, so including or excluding the tip of the cone 
(contrast the products in eq. (14. 2 p and eq. (14. 3p ) is immaterial. In this case the condition 
(I4.34P may equally well be taken to specify the exact number of traces or the maximal number 
of traces. We call this regime the classical phase, to distinguish it from the Bose-Einstein 
phase below. For completeness, we quote the usual formula for the entropy in the classical 
phase: 

S{n, M) = log Z + I3{n, Vf{n)) - aM . (4.43) 

4.2.2 Bose-Einstein phase with trivial point included 

When M > Mcrnin), a macroscopic proportion of trace factors condense on a small set of 
points and the system enters the Bose-Einstein phase. Depending on whether one works 
with partition function (14. 2 p or (14.30 . that set of points is either the tip of the cone or the 
points nearest to it. 

We begin by including the trivial point and consider partition function (14.21) . This cor- 
responds to fixing the maximal number of traces at M. The requisite counting is given 
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by the generating function QM^i) in [T]. In this case, the condensate at corresponds to 
M — Merit (^^) > excess trace factors, which the system dumps in the immaterial ground 
state. The physical meaning of the Bose-Einstein condensate is clear: it is the phase where 
it is more entropically favorable for the system to decrease the number of (non-trivial) trace 
factors. 

Following the standard treatment of Bose-Einstein condensates, we have 



B{n)T{d+ l)C{d+l) 



in, V fin)) = (4.44) 



where in the second line we assumed Mbc, the population of the condensate, to be large. 
The entropy takes the form: 



S{n,M) = log Z + /5(r2, V/(n)) -a(M-M, 



be J 



\ogZ + /5(n, Vfin)) + + 0{M^^') . (4.46) 

IVl] 



be 



In the large Mbc limit, only the first two terms remain. These, of course, correspond to the 
entropy contained in the MCTit(?^) excited particles (non-trivial traces.) Therefore, f l4.46p is 
expected to agree with (13.201) by the standard arguments establishing the equivalence of the 
microcanonical and grand canonical analyzes in the thermodynamic limit. 

4.2.3 Bose-Einstein phase with trivial point excluded 

We now work with the partition function 

mec*\o mec*\o 

This case corresponds to fixing the exact number of traces at M > Mcritin). 
Each state in the system is characterized by a choice of such that 

m = n . (4.48) 

As seen from eqs. (14.101 14.351) . the occupations controlled by the quantity y{fh) = 

(m, V/(fi)), which plays the role of the energy of the system. With the tip of the cone 
excluded, the Bose-Einstein condensate will accumulate on those integral points in C* for 
which y{rn) is the smallest. We denote this minimal value with yQ. The structure of a state 
as a vector partition (14.481) descends to a partition of the quantity y{n) = (n, V/(n)): 

yin) = ^ y{m) = Mbc 2/0 + ^ '^m y{m) . (4.49) 
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In the last equality we have grouped the contributions to y{n) of the Bose-Einstein condensate 
and of the excited particles, respectively. 

Let < c < 1 be the proportion of yln) contributed by the second term in (14.491) : 

^^n^?/(m) = cy{n) ^ ''^^ rirn rh = c n (4.50) 

m^bc m^bc 

Mbc2/o = (l-c)y(n) (4.51) 

The condensate does not begin to set in until the population of all the excited particles 
reaches the critical value. The latter is given by 

M„it(cn) = c^M„it(n) (4.52) 

since the charge vectors of the excited particles form a partition of cn as in eq. (I4.50p . The 
scaling of Merit (cn) with c follows from the definition (14.411) and the fact that f{n) is a 
(i-dimensional volume. Then imposing the condition that the total number of traces be M 
is equivalent to requiring: 

M = (1 - c)^ + c^M„it(n) . (4.53) 

yo 

Using the fact that eq. (I4.50p is a partition of cn with a critical number of traces, one may 
feed the solution of fl4.53p into eq. fl4.20p to get: 

^ B{n)T{d+l)ad+l) 
(n, V/(n)) 

As before, the chemical potential in the large Mbc approximation is given by 

«^-lf^ = -(l-c)-^^. (4.55) 
Mbc y{n) 

This yields for the entropy 

S{n,M) = logZ + c/3(n, V/(n)) -a(M-Mbc) 

^ Bin)r{d+l)ad+l) yoMcHt(r^) 

1 - c y{n) ■ ^ ' ' 

At fixed n, the limit of large M is the limit of small c and consequently, from eq. (I4.54p . 
large (3. In the latter case, we see that every term in (I4.56P goes to 0. This, of course, agrees 
with the intuition from the microcanonical counting: when the number of summands grows 
too large, the number of partitions of n quickly decreases, reaching when M oc \n\. 
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4.3 Typical configurations of dual giant gravitons 

As in Sec. [3l the results of the present section afford an interpretation in terms of dual 
giant graviton configurations. Indeed, the partition function (14.21) has been written down 
in [4] specifically to count configurations of dual giant gravitons. In particular, given the 
total momentum P on the maximal torus of X, a typical distribution of momenta of the 
individual dual giant gravitons is recovered using eq. (14.101) and its generalization as in eq. 
(I4.35p . This is done, as in Sec. [3l by selecting a special cycle = o"v/(p) ^ ^- then 
sets the momentum component along o'yji^p) for each dual giant graviton using the (grand) 
canonical ensemble. The distribution of 3-momenta of the individual dual giant gravitons 
follows from this procedure for essentially statistical reasons: the cycle c'"v/(p) selected in 
such a way that momentum along it splits into the three fundamental cycles of the torus 
precisely in the ratio given by P. 

The different ensembles considered above have natural interpretations in the language of 
dual giant gravitons. For completeness, we dilate on this below: 

1. The analysis of Sec. 14. II applies to the case where one sets the total momentum of the 
dual giant gravitons to P ~ fi but leaves the total number of dual giant gravitons 
unrestricted. The ensemble is then left to select its own, optimal number of dual giant 
gravitons. 

2. The analysis of Sec. 14.21 applies when one wishes to fix the total momentum of the dual 
giant gravitons to P ~ n and the total number of dual giant gravitons to M. The 
previous expressions are modified via eq. (I4.35p . Several subcases are distinguished 
depending on the relation between M and Merit (P), eq. (14.411) : 

2.1 If M < Merit (P) then eq. (14.341) sets the parameter a in a simple manner. In this 
case M can be thought of as fixing the exact number of dual giant gravitons, or 
the maximal number of dual giant gravitons. 

2.2 If M > Mcrit(P) fixes the maximal number of dual giant gravitons, the results of 
Sec. 14.2.21 apply. In that case, the number of dual giants that maximizes entropy 
equals Mcrit(P). 

2.3 If M > Afcrit(-P) fixes the exact number of dual giant gravitons, the results of 
Sec. E3] apply. 

4.4 Summary: The universal features of heavy operators 

In all the theories we have considered, multi-trace operators and their single-trace factors are 
charged under three f/(l)s. We have shown that for a given total charge n, the distribution 
of charges of single-trace factors under a certain U{1) is universal. This U{1) is selected 
by a function f{n) on the dual toric cone which measures the minimal volume in the cone 
cut out by a plane through n. Given a triplet of charges rh, the charge of this particular 
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f/(l) is computed by the inner product t(m) = (m, V f{n)). Then the fraction of single trace 
factors carrying U(l) charge Ci |n|^/^'^^"'^^ < t(m) < C2 \n\^^^^~^^^ for any constants Ci and C2 
will be universal. The distribution of charges under the two residual U{l)s is not universal. 
In addition the total number of traces in a heavy multi-trace operator is proportional to 
rjij^g proportionality factor depends only on the dimension d of the toric variety. 
One concern about this result is that generic multi-trace operators may not be orthog- 
onal and hence we may wonder whether an orthogonal basis enjoys the above universality 
properties. To examine this, we can imagine starting with one of the above operators and 
then systematically producing an orthogonal basis by a procedure like Gram-Schmidt. By 
the arguments given above, almost all the summands appearing in the orthogonal linear 
combinations will have the same trace structure. In this way, it appears that the universal- 
ity will extend to an orthogonal basis. In any case, the dual giant graviton interpretation 
is free of this caveat. In that case, the universal properties include: (a) the number of dual 
giants in a typical state with total momentum P, and (b) the fraction of the giants carrying 
a fixed amount of momentum around a particular cycle denoted Cyfi-pv 



5 Dual black hole geometries 

In Sec. [3]we have enumerated heavy scalar mesonic BPS operators with a given charge vector 
n. The structure of such operators exhibits many universal features, which were described 
in Sec. IH either in the trace basis or using the dual giant graviton interpretation. The 
AdS/CFT correspondence anticipates that the dual description of such states is in terms of 
black hole geometries asymptotic to AdSs. We now turn to reinterpret our previous results 
in terms of black hole physics. 



5.1 Dual black holes have vanishing horizons 

Here we argue that the dual black hole geometries have semi-classically vanishing horizons. 
Since by the Bekenstein relation the area is proportional to the entropy, S = A/{4Gn), in 
our case we have 

Vol(X) pI 

^grav. T^TF' (5-1) 

> ''p 

where ph is the radius of the horizon and Vol(X) is the volume of the transverse Sasaki- 
Einstein manifold X. This volume has been computed in [15], and is proportional to the 
volume of the 5-sphere as 

7777^ = = h{ai), 5.2 

where the constant of proportionality, h{ai), is the ratio of the central charges of the two 
dual SCFTs. It was shown in [9] that the volume of X is proportional to the volume of the 
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polytope A in the dual toric cone C* bounded by the plane 2(6, n) = 1: 

Vol(X) = 6(27r)3 Vol(A) = 2(27r)V((0, 0, ^)) . (5.3) 

The second equality follows from the discussion above eq. fl3.2ip . As /55((0,0, ^)) = ^, we 
have: 

%,) = 48/((0,0,-)). (5.4) 

Using the relation between string length and Planck length, Ip/lg ~ gl^^, and the volume 
of the 5-sphere, Vol(S'^) ~ {gsN)^^^l^, we solve for the horizon radius in string units as 

I J ^[ij VoliX) [ij ^''^''■^ ig,N)ihia,)^ h{a,) ' ^'''^ 

This indicates that the radius of the horizon is nonzero in string units only if the entropy of 
the black hole scales as A^^ or faster, where N is the number of D-branes giving rise to the 
SCFT. 

Naively, one would expect the dual geometries to correspond to black hole microstates 
once the conformal dimension is of the order A^^. This is because in other contexts black holes 
have been described in terms of a collection of A^ D-branes, each with conformal dimension 
of the order A^ [16]. These states do not have macroscopic horizons, however. To see this, 
consider any lattice site n in the dual toric cone whose conformal dimension is A. The 
volume of any pyramid cut out of the dual toric cone with a plane through n will grow as 
|n|^. Therefore f{n) oc N^'^ and by eq. IKM : 

S{n) oc A^^ = N^-^ < A^^ _ (^5_g) 

We conclude that the entropy never grows sufficiently fast to produce a horizon. 

One might object that the above argument starts out treating A^ as a large but finite pa- 
rameter while the calculations leading to eq. fl3.20p assumed infinite A^. The latter statement 
is an immediate consequence of working with dual toric cones, which extend along infinite 
rays. By contrast, in theories defined on finite stacks of D-branes, the ring of BPS operators 
will be subject to some non-trivial relations, due to the fact that for a rank A^ matrix, the 
operator Tr {X^~^^) can be written in terms of lower dimensional traces Tr(X), . . ., Tr(X^). 
These relations are generated by Newton's identities applied to minimal loop operators and 
begin to kick in at distances of order A^ away from the tip of the cone. As a result, lattice 
points with |f2| > A^ in general fail to produce algebraically independent trace factors. We 
conclude that a finite A^ may be thought of as defining a truncation of the cone, so that all 
trace factors (partition summands) must fall within a rough distance A^ of the tip. 

Happily, below eq. (I4.32p we argued that, in almost all partitions, almost all summands 
satisfy \rn\ oc jnl'^+i oc N ^+1 , so that |m| -C A^ for (i > 1. Therefore, a generic partition does 
not see the truncation by finite A^. This establishes the consistency of the argument leading 
to (15.61) . As an extra observation, notice that the inequality (15.60 is improved as d increases. 
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Finite horizons: Equation fl5.5p tells us how the entropy of a black hole solution should 
scale for the black hole to have a macroscopic horizon. If we choose our conformal dimension 
to scale as A ~ A^''', we see that for a macroscopic horizon we need 



^ 5 ~ ~ AsTT ~ A^dTT ^ ^ = '^'^^ + ^) ^ !i for = 3. (5.7) 

a 3 



Thus the states of interest to us have conformal dimensions A ~ N^^'^. Luckily, this is 
still not large enough to invalidate the argument above: in almost all partitions almost all 
summands satisfy |m| ~ ~ N~ ~ A^^^"^, so that \rfi\ <^ N. It is interesting to note 

that 7 decreases as d increases in keeping with the fact that less supersymmetry makes it 
easier to develop a horizon!^ 

Gravity duals of heavy operators: For states carrying non-vanishing R-charge in the 
superconformal field theory, the dual metrics were given in [T7|, who used the consistent 
truncation discovered in [18]. This R-charge is realized in supergravity by turning on angular 
momentum along the cycle defined by the Reeb vector. Denoting by the Killing vector 
field along the Reeb vector, the metric reads 

ds"^ = -^H-^Fdt^ + H[F-^dr^ + r^dsl^] + ds^ + + Af, (5.8) 
where F is a four dimensional Kiihler-Einstein manifold, and 



H = 1^ 



F = 1 + (5.9) 

A = -H^Ht. 
2 

The gauge field A is responsible for the non-vanishing R-charge. These metrics have mani- 
festly vanishing horizons, in agreement with our field theory counting and stretched horizon 
considerations. They carry only R-charge and are not charged under any extra U(l)s. This 
is equivalent to choosing our charge vector n to point in the z-direction (in the standard 
presentation), as dictated by fl3.23p . 

For 7W = C'^, (singular) supersymmetric configurations with vanishing horizons have 
been known to exist for any value of the three U(l) R-charges, as long as the system remains 
BPS [19]. This statement matches our observations for any charge vector n. 

Since our methods and results apply to general toric varieties M. with arbitrary U(l) 
charges turned on, any comparison with gravity would require the addition of these extra 
quantum numbers (nonzero fiavor charges). To our knowledge such solutions have not been 
written down in the literature. We claim that the relevant geometries, when discovered, will 
be characterized by vanishing horizons. 



It is unclear to us whether such heavy operators will have a classical spacetime description in terms of 
an asymptotically AdS^ x 5*^ metric, or whether this asymptotics will be modified. 
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6 Discussion 



In this paper we developed a geometric way of counting scalar mesonic operators in A/" = 1 
toric quiver gauge theories based on the one-to-one correspondence between (equivalence 
classes of) gauge invariant operators and lattice points in the corresponding dual toric cone. 
This correspondence also allowed the recent development of the so-called plethystic program 
for counting the same quantities [1]. 

Our geometric considerations allowed us to analyze the structure of mesonic operators in 
the large charge limit using statistical arguments. This was accomplished by working with 
the appropriate partition functions (ensembles), and realizing that the set of states having 
the same occupation numbers form planar surfaces in the dual toric cone, inducing a natural 
grading on it. This connection allowed us to find the direction in the dual toric cone that 
had the most rapid growth of entropy. 

In the large charge limit, the grading becomes a continuous variable y. We were able to 
use this to encode the structure of gauge invariant operators in terms of curves Xd{y), whose 
properties only depend on the dimension d of the initial toric variety!^ Here, Xd{y) stands 
for the ratio of traces in a given multi-trace operators having grading value y < y. Similar 
mathematical functions control the distribution of all f/(l) charges among the traces of the 
heavy multi-trace operator. 

Our statistical analysis shows that there exists a universal structure, shared by almost 
all heavy operators. In this sense, given an arbitrary constraint like fixed R-charge, there 
exists a typical state close to which almost all other states lie; these can be thought of as 
master states in the field theory. This provides a generalization, for any M = 1 toric quiver 
gauge theory, of the results and philosophy reported in [13] for the half-BPS sector of A/" = 4 



We mainly worked on the field theory side, and according to the AdS/CFT duality 
appropriate combinations of these operators are expected to be dual to microstate geometries 
on the gravity side. Unfortunately, the exact form of these geometries is not known for a 
general quiver0 Should these geometries be discovered in the future, our methods could be 
used to analyze the emergence of semi-classical gravity from the fundamental microstates in 



Even if we would know the explicit candidate microstate geometries, we know that their 
coarse-grained description would be interpreted in terms of distributions of giant gravitons or 
dual giants in AdSs x X. All our results can actually be reformulated in terms of dual giants 

In other words, the details about the precise nature of the toric variety are suppressed by negative 
powers of the charge, or equivalently, are suppressed in N . 

Recently, for the special case where the corresponding toric variety is C^, the equations describing the 
moduli space of such candidate microstate geometries were constructed in [20] . 

Some recent work in the literature concerns the emergence of locality and spacetime from a study of 
the strong coupling dynamics of the gauge theory dual to the conifold [12] . It would be of interest to apply 
these techniques both to our heavy operators and to more general N = 1 setups. 



SYM. 




38 



due to the isomorphism of Hilbert spaces discussed in Section I2.4i^1 Since our universahty 
conclusions also apply in the description in terms of dual giant gravitons, this strongly 
suggests that the physical picture of black hole microstates advocated in [13] will also hold 
for all the systems described in this paper. 

We have only worked at a classical level, and giants and dual giants are quantum me- 
chanical objects. Thus, to establish the validity of our analysis, we need to show that the 
wavefunction of the configuration of D-branes is strongly localized, indicating that the sys- 
tem has a classical description. We established this localization for the case of giant gravitons 
on in appendix [B] and expect it to hold in a general setting. 

In Section [5] we showed that to produce a macroscopic horizon, we would have to consider 
BPS mesonic states with conformal dimensions of at least of order N^^^, for the case of a 
three dimensional toric variety, i.e., for string theory on AdSs x X. It would be interesting 
to understand what these states look like on the gravity side, and why such large conformal 
dimensions are needed to create horizons. Naively one would expect a black hole to arise from 
a collection of 0{N) D-branes and thus have a conformal dimension N ■ N N'^ <C iV®/^. 
However, whatever the gravitational description, the structure of these states can be studied 
using our methods, barring a breakdown of the canonical ensemble it self 1^ 

We would also like to point out that the case of a four (complex) dimensional toric variety 
might be interesting to study. In this case, the conformal dimension required to produce 
a macroscopic horizon is smaller, N^^"^, and the field theory is simpler. This is because 
increasing the (real) dimension of the toric variety by two, correspondingly decreases the 
dimension of the transverse Minkowski space by two, and after the near horizon limit this 
yields string theory on AdSa x X, where X is is the seven- dimensional base of the toric 
variety. This is naturally expected to be dual to a (1 + l)-dimensional CFT. 

Finally, we argued that even though the trace basis is not orthogonal for large values 
of the U(l) charges, the typicality we found would also be present in an orthogonal basis. 
However, it would be interesting to work explicitly with such an orthogonal basis, analogous 
to the basis of Schur polynomials for the half-BPS sector of A/" = 4 super- Yang-Mills |23j . 
Such a basis was recently provided in [21] for the chiral ring of A/" = 4 super- Yang-Mills. 
It would be interesting to extend that analysis and techniques to chiral rings of arbitrary 
A/" = 1 super conformal field theories. 

We expect that last point to be relevant for a comparison with gravity. This is because 
a generic heavy multi-trace chiral operator has no natural giant graviton interpretation [21] 
and we already know that Schur polynomials played a crucial role in matching the analysis 
of typicality in gauge theory with gravity [T3] . 

As discussed in appendix[Cl the existence of the identity between partition functions (|2.13p suggests that 
a similar interpretation in terms of giants should also be possible, and would be more naturally formulated 
in terms of fermions. 

The ensemble might conceivably break down due to a Hagedorn-like growth in the number of states; it 
might be interesting to study this possibility further. 
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A Exact multiplicities of operators with irrational R- 
charges 

In Section 13.11 we showed that the combined entropy of operators with co-planar charge 
vectors n is a function (eq. (13. lip ) of the volume which the plane cuts out of the dual 
toric cone. The simplest application of this statement is to evaluate the total entropy of 
all operators of fixed R-charge Rq. The charge vectors satisfying this constraint live on the 
plane 

{b,n) = ^Ro. (A.l) 

The volume of C* bounded by the plane determines the total entropy via eq. (13. lip . 
The above statement is valid in the following two cases: 

• When the Reeb vector b contains only one Z-linearly independent component, 

• When one wishes to impose flA.ip up to some small finite accuracy e. 

Suppose, however, that among the components of the Reeb vector there are k > 1 Z-linearly 
independent ones. Denote them by rj, z = l,...,k. One example of that is dPi, where 
b = (0, 4 — vT3, 3) and 4 — vT3, 3 are Z-linearly independent. Suppose further that one 
is interested in constraining the R-charge at exactly Rq. We may decompose Rq and b into 
summands proportional to rj: 

b = niblblbl) + ... + rkiblblb',) 

Ro = nRl + ... + ruR^o- (A.2) 

If no decomposition of i?o is possible, then no state with R = Rq exists and computing 
the entropy is meaningless. Otherwise, the constraint (lA.ip decomposes into k independent 
ones: 

{b\n) = -R, z = l,...,A;. (A.3) 

Therefore the locus of charge vectors exactly matching Rq is given by the region of intersection 
of C* and k planes (lA.3p . In other words, it is the section of a certain [d — A;)-hyperplane 



contained in the dual toric cone. 

For (i = 3, three cases are distinguished: 
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1. k = 1. The locus R = Rq is a plane. The combined entropy is determined by the 
volume of C* contained under the plane as explained in Section 13.11 

2. k = 3. The locus R = Rq is a point. The entropy is determined by the minimal 
volume of C* that may be cut out with one plane through the point R = Rq. This was 
explained in Section [3121 

3. k = 2. The locus R = Rq is a. line. This case, along with all others with 1 < k < d, is 
treated in the next subsection. 

A.l Multiplicities of collinear vectors 

The argument of Section [3l2] is readily adapted to 1 < k < d. The locus of points of interest is 
a (c? — /c)-hyperplane satisfying k equations (]A.3p . Consider the family of {d— l)-hyperplanes 

containing the locus R = Rq. Here s is the normal to the hyperplane, which is specified 
by A; — 1 free parameters. In particular, s is a linear combination of the vectors 6* whose 
normalization is immaterial and may be scaled to 1. 

For each such s, consider the family of {d — l)-planes parallel to it and proceed as in 
Section [3TT1 The total entropy of the points {R = Rq} is bounded from above by 

S{Rq) < min S{p,) . (A.4) 

Here S{pg) denotes the total multiplicity of all charge vectors contained in pg, which may be 
obtained from eq. (13.111) . Assume that every {d — /c)-hyperplane {R = Rq} in p^^i„(_Ro) has 
its own minimal Smin(-Ro)- Since the region {R = Rq} is contained in p^„,i„(_Ro)5 its minimum 
must be strictly smaller than the entropy of the {d — l)-plane p?„,i„{_Ro)- This applies to every 
locus {R = Rq ^ Rq}, which implies that: 

SiRo) ^ min Sipg) . (A.5) 

This equation applies to all 1 < A; < and naturally generalizes eq. fl3.13p . 

As k grows, one is able to derive more information from a single exact measurement of the 
R-charge. Such successive specifications narrow down the microcanonical ensemble of states 
{R = Rq}, which should result in a decrease of the multiplicities computed in eq. (1A.5I) . We 
may see that this is indeed the case from the schematic relationship 

S{RQ)\k=d ~ min S{pg) 
< 5(i?o)|fc=d-i ~ min S{pg) < ... 

<S{Rq)\,=, = Sip^). (A.6) 

The multiplicities obtained from specifying the charge vector n with a growing number of 
constraints differ only in the domain of minimization in eq. ( 1A.5I) . In particular, they are of 
the same order of magnitude. 
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B Mapping to the typical giant configuration 



In Section 14.31 we interpreted our analysis as giving a typical configuration of dual giant 
gravitons. However, all these configurations are at heart quantum mechanical, and to analyze 
typicality we should verify that the wavefunction of the typical configuration is sharply 
peaked, thereby allowing one to associate it to a semiclassical collection of D-branes. In this 
appendix we will argue that this indeed is the case: the typical states constructed above can 
be associated to a typical configuration of giant gravitons on the gravity side. For simplicity, 
we will deal with giant gravitons on S^; we expect the conclusions to be valid in a more 
general setting as well. 

It was shown by Mikhailov [25] and Beasley [11] that these brane configurations can be 
given in terms of zero-loci of polynomials of the form 



where n runs over N^. The worldvolumes of the branes are given by the intersection of the 
surface P{z) =0 with the unit in C^, where the is identified with the one in AdSs x S^. 
The phase space of brane configurations is thus spanned by {c„^„2„3} and is topologically 
CP°°. 

This space was quantized in [26] using geometric quantization. First one defined a new 
set of coordinates that eliminated the pathological regions in coordinates {c^^^jns}, i-e., 
eliminated polynomials that don't intersect the unit and identified polynomials that have 
the same intersection with the S^. The new coordinate system is given as 

where f\ni+n2+n3\ is a smooth function. Upon quantization the Hilbert space of brane config- 
urations was found to be the space of degree N holomorphic polynomials in Wnin2n3 , which is 
isomorphic to the Hilbert space of N free bosons in a three dimensional harmonic oscillator, 
and therefore has the partition function 

Z = Y\ r- (B.3) 

Above, we've restricted the highest excitation any particle can have to be given by the 
constant N^, physically this corresponds to fixing the the number of giant gravitons |13j . 
This is because the highest excitation gives the degree of the polynomial P{zi, Z2, z^) and 
therefore the maximal number of D-branes. In the following we'll also work in the high 
temperature regime, (3 —>■ 0. 

We can now write down the wave function of the typical state. Since the energy levels 
factorize, we have 

<ij'yp''^\w) = Y[^j{w), (B.4) 

j=0 
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where ipj is the contribution from particles with energy j to the total wavefunction. Since 
the wave function has to be symmetric both under the exchange of particles and under the 
exchange of dimensions, we can write 

j 



E 



nx,ny,n2=0 



Since the moduli space is CP°°, the coordinates have to satisfy 



x,ny,nz \'^^x,ny,nz 



(B.5) 



1. 



If we define the contribution from each energy level j as Rj, we can write this condition as 



nx,ny,nz=0 
nx+ny+nz=j 



R], with ^R] = l. 

j=0 



For any given Rj, we know that the maximum of | Yl''n 



.^,ny,n2=0 '^nx,ny,n, 
nx+ny+nz=j 



(B.6) 



occurs when all 



the w's are equal, and since the number of states at energy level j is ("'^^), we see that the 
correct value is 

y ni + n2 + n3=j. (B.7) 



(Up to an overall phase that is not relevant to us.) Thus at the maximum the wavefunction 
can be written as 



A^c r 
j=0 



J + 2 
2 



(B.8) 



i=o 



where the normalization constant C is not important and will be dropped. We wish to 
maximize |\E'p with the constraint YlR] — ^^"^ to do this we eliminate Rq by solving the 
constraint, yielding 



2ki 

i ■ 



k=l 



1=1 



After some algebra one finds the maximum of this at 



Rj 




iV' 



(B.9) 



:b.io) 



which by construction satisfies R^ ~ ^- This is the only local extremum, and since 
the wavefunction vanishes on the 'edges' where wj = for some j, this is also the global 
maximum. 

Finally, we compute the second derivative of the wavefunction at the maximum. After 
some algebra this yields 



1 d^l^iwW 



|^(w)|2 dRjdRi 



(B.ll) 
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Since in the large temperature limit /5 — > the free energy is minimized when the entropy is 
maximized, we see that for the typical state the particles are evenly distributed among all 
possible state£§ , and therefore the number of particles on energy level j is kj oc (•'^^) • Thus 
we see that in the limit — > oo the wavefunction is sharply peaked around the maximum, 
and we can associate the typical state with the location of the maximum in the phase space, 
^max^ The map flB.2p then guarantees that there exists a corresponding point c^""^ specifying 
the typical configuration of giant gravitons. We expect this configuration to correspond to 
the superstar [16], but have been unable to verify this as the explicit form of the funtion 
/|ni+n2+n3|(c) in (lB.2p is not known. 



C Limiting Curves and Amoebae? 

There have recently been developments in topological strings using the statistical physics 
of so-called melting crystals [271 [28] . Analyses of typical states in black hole physics have 
involved the use of similar statistical mechanics techniques. In this brief appendix we would 
like to comment on the similarities and differences between the two and make some specu- 
lative remarks. 

In the melting crystal studies one considers geometric quantization of a single particle 
in C"*. Its spectrum is that of a harmonic oscillator in m dimensions, and as such, the 
associated partition function is related to a classical partition function. For example, for 
C'^, there is a one-to-one map between quantum harmonic oscillator states and 3-D Young 
tableaux. If we associate each tableaux with a distribution of boxes in an octant of M^, 
the statistical mechanics provides the limiting shape for the boundary of the crystal as a 
function of the string coupling, which plays the role of the temperature in this set-up. 

This discussion is conceptually very reminiscent of the derivation of the typical state in 
the half-BPS sector of A/" = 4 SYM. There is a one-to-one correspondence between N -particle 
fermionic quantum states in a 1-dimensional harmonic oscillator and 2-dimensional Young 
tableaux. Statistical mechanics and a large limit allow the derivation of a limit shape 
y{x) which describes the average excitation {y) of the particle x. 

Though the techniques are the same, the physics looks rather different. The plethystic 
methods and the geometric quantization of the classical moduli space of dual giants provide a 
clue as to where a relation may exist. Let us consider the isomorphism, which is a cornerstone 
of the plethystic program: 

gN{t,M)=gi{t,M''/SN). (C.l) 

This states that the number of mesonic scalar chiral primary gauge invariant multi-trace 
operators constructed out of building blocks that define M. is equal to the number of single 
trace operators in a different manifold, the symmetric product M.^ /Sn. In terms of dual 
giant gravitons, this has a very natural interpretation: the Hilbert space of dual giants is 
equivalent to the symmetric space of the Hilbert space spanned by a single giant, which is 

This is also easy to derive from (|B.3I) more rigorously. 
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the set of holomorphic functions in Ai. Therefore, fIC.ll) provides a bridge between partition 
functions of particles and partition functions of a single particle in a different space. 

Consider the particular case oi A4 = C The Hilbert space of a single particle is that 
of a harmonic oscillator in one dimension. When we consider the Hilbert space of a single 
particle in /Sn, all quantum numbers in each of the copies appearing in the symmetric 
product have to be different: this is effectively implementing the Pauli exclusion principle, 
and so we conclude that we are equivalently dealing with fermions in a one dimensional 
harmonic oscillator. The limiting curve, exhibited in eq. (44) of [13], is given by 

ae^^ + be-f^y = 1 , a = e"^^, b=l-a (C.2) 

and reproduces the curve xi{y) of eq. fl4.29p up to a linear re-definition of coordinates. 
Inspecting Section 4.1 of [28] we see the emergence of flC.2|) therein. Equivalently, in the 
language of [29], ( ]C.2I) is the boundary of the amoeba of C"^, which, in turn, is equivalent to 



the limiting shape of constructing three-dimensional Young tableaux. 

We can now understand the space where the limiting curve y{x) lives. Indeed, x can be 
understood as parameterizing the copy x in the symmetric product C^/ (or the particle 
X in the dual giant language) whereas y is just describing the excitation of that copy (par- 
ticle). In the semiclassical limit x is a continuum variable, and we obtain a curve on a two 
dimensional plane. This curve describes a single particle living not in C^, but rather in the 
symmetric product C^/^at with N oo. 

To recapitulate, we have that (1) the curve of typicality of half-BPS states of A/" = 4 SYM 
identifies with (2) the limiting curve of the topological A-model on C'^. We point out that 
the underlying geometries of the two differ: the fundamental generating function (i.e., for the 
single-trace operator) of (1), in the language of the plethystic program, is the Hilbert series 
for C while the curve in (2) is the amoeba for C'^. Nevertheless, the emergence of the same 
curve in two different counting problems suggests that there may be a deeper connection, 
possibly persisting to other geometries, and hence other A/" = 1 theories, which we should 
explore. 
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